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ABSTRACT 

Inflationary models predict a definite, model independent, angular dependence 
for the three-point correlation function of AT/T at large angles (> 1°) which we 
calculate. The overall amplitude is model dependent and generically unobservably 
small, but may be large in some specific models. We compare our results with 
other models of nongaussian fluctuations. 

Subject headings: cosmology: theory — large-scale structure of universe 



Recent results from COBE (Smoot et al. 1992) have provided strong evidence for 
fluctuations AT in the cosmic microwave background radiation (CMBR) on large angular 
scales 1°). The two-point correlation function of AT is consistent with a Harrison- 
Zel'dovich spectrum, as predicted by inflationary cosmologies (Guth 1981). Here we point 
out that inflation results in a deflnite and model independent prediction for the angular 
dependence of the three-point function as well. The overall amplitude of the three-point 
function, however, is model dependent. It is possible to construct models (discussed below) 
where this amplitude would be large enough to be seen in the COBE data, but generically 
it is expected to be too small to be observable. 

For large angular scales, AT is predicted to be (Peebles 1982; Abbott & Wise 1984; 
Bond & Efstathiou 1987) 
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AT(a)=^^J Qjd\me''-\ (1) 

where a is a unit vector on the sky, Q is the predicted RMS quadrupole amplitude, and 
the two-point correlation function of the random variable ^(k) is given by 

{m)i{k2))^kih\ki + k2), (2) 

where 5^(k) is the three-dimensional Dirac delta function. To study the two-point corre- 
lation function of AT in a realistic setting, we must flrst expand AT in Legendre polyno- 
mials, remove the monopole, dipole, and (following Smoot et al. 1992) quadrupole terms, 
and weight the remaining terms to account for the flnite beam width a of the antennas. 
This results in (Bond & Efstathiou 1987) 

00 

^ik-a _^ J2i\2l + 1) Wpiik) m-a) , (3) 
/=3 

where k = \k\^ k = k/k^ Jlik) is a spherical Bessel function, and Wi ^ exp[— ^(/ + ^)^a"^]; 
Smoot et al. (1992) take a — 3.2°. This implies that the two-point correlation function for 
AT takes the form 

(AT(a)AT(6)> = |g2^2(a,S), (4) 

where 

oc 

C2(a,6) = J2i'^l + l)CiW^Piia-b) , (5) 

with Ci = 2 dkjf{k)/k = + 1). For reference, €2(0,, b) is plotted as a function of 
the separation angle a = cos~^{d-b) in Figure (1). 



In inflationary models, the random variable ^{k) has its origin in the fluctuations 
of a quantum field (p. Working in momentum space, and assuming (p is massless and 
noninteracting, during the de Sitter epoch the two-point correlation function of ip is 

(^(^1, t)^(^2, t)> = A7r^H^k^'-^5'\ki + k2) (6) 

in the long wavelength limit (Bunch & Davies 1978). Here r is conformal time and H is 
the Hubble parameter, related to the (constant) scalar potential Vq via = 87rVo/3. (We 
take h = c = G = 1.) The three-point correlation function vanishes identically. However, 
if the field (p has any sort of interaction (with itself or with other fields) , nontrivial higher- 
point correlation functions will appear, and these will be passed on to the temperature 
perturbations (Allen, Grinstein, & Wise 1987). The simplest possible interaction, and one 
which will generically be present, is a modification of the scalar potential from a constant 
to V{(p) = VQ + ^fi(fi^. To leading order in this interaction implies a three-point function 
which can be calculated via standard methods of quantum field theory in curved space 
(Birrell & Davies 1982). In the long wavelength limit, we find 

{(fi{ki, T)(fi{k2, rMh, t)> = §71^/^2 [kik2k3)-''^ F{ki, k2, ks) 

xS^{ki + k2 + h) , (7) 

where 

F{ki, k2, ks) = (kl + kl + kl) [log(A;T|T|) + 7] 

- {kj + k^ + kl)k^ + kik2k2, . (8) 

Here ki = \ki\ and kT = ki + k2 + k^: 7 is Euler's constant. The factor of log(/cT|T|) 
in equation (8) is minus the number of e-foldings from the time at which fluctuations of 
wavenumber /ct first passed outside the horizon until the end of inflation. For the relevant 
length scales, one typically finds log(A;T|T|) — —60 (e.g., Bardeen, Steinhardt, & Turner 
1982), and so the term with this factor dominates over the others. From here on we will 
approximate F{ki, k2, k^) as 

Fiki,k2,k3)c^-Pikl + kl + kl) (9) 

and treat /3 ~ 60 as independent of k^. A similar logarithmic dependence on k^ (with a 
different origin) is also present in equation (2), and is also routinely ignored. 

The fiuctuations in (p become fiuctuations in the CMBR in one of two ways. If the 
energy stored in (p is eventually converted to ordinary matter and energy (including the 
CMBR), then the fluctuations are adiabatic, and the correlation functions of ^{k) are the 
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same (up to a sign) as those of (p{k)/{A7r^H^)^/-^. In this case we typically find /i/H ^ 
10~^ (Hodges et al. 1991), but larger values may be possible if we fine-tune the scalar 
potential. If, instead, the energy stored in (p is eventually converted to dark matter, then 
the fiuctuations are isocurvature (Linde 1984), and the correlation functions of ^{k) are 
again the same (up to a sign) as those of (p{k) / {Att^ H'^)^/'^ (Allen et al. 1987). In this case 
there are no obvious restrictions on fi/H. 

Neglecting, for the moment, the effects of finite beam width and the need to subtract 

the Z = 0, 1, 2 terms in AT, we can combine Eqs. (1), (2), (6), (7), and (9) to compute an 

"unsubtracted" three-point correlation function for AT. We find 

-|/o\3/2 ^ 
(AT(a)AT(S)AT(c)>^^3^b = ^ - [^) P^Q' Cs{a, b, c) , (10) 

where £ = ±1 is a model dependent sign, and 

C^{a, b, c) = ^ ^ I d\i d% d% S'^{ki + k2 + h) 



3(27r)^ 



X (fcf 3^2-3 + . . .) e^ikv^+h-b+h-c) 



3(27r)2 

C2(a, c)C'2(6, c) + C2{b, a)C2{c, a) + ^2(6, 6)^2(0, b) 



(11) 



In the last line, we have defined an "unsubtracted" dimensionless two-point correlation 
function 



C2(a,6) = 7^ / d\k-^e'^<^-^^ 
2tt J 



= ^{2l + l)CiPi{d-b). (12) 

1=0 

Equation (11) is a beautifully simple formula, but clearly is only a formal relation, since the 
monopole term in C2 is infinite. Equation (11) must be corrected to remove the / = 0, 1, 2 
terms, and to account for finite beam width. This is straightforward, and ultimately yields 

(AT(a) AT(6)AT(c)> = ^ - (^5 j P^Q^ C^ia, 6, c) , (13) 

where 

00 

C3{a,b,c) = l J2 i2j + l)i2k + l)i2l + l)iCjCk + CkCi + CiCj) 
j,k,l=3 

xWjWkWifjkiia,b,c), (14) 
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and we have defined 



1 



/ 



(15) 



Here dVtn denotes integration over the unit vector n. If desired, f jkiid, 6, c) can be expanded 
in spherical harmonics of each of the three unit vectors by making use of standard relations 
[for example, eqs. (9.79) and (16.90) of Merzbacher (1970)], but the resulting formula for 
the coefficients is unwieldy and will not be presented here. 

C^{djb,c) is plotted in Figures (2) and (3) for an equilateral triangle on the sky: 
d-b = b-c = c-d = cos a. Surprisingly, for this case it is given to a good approximation by 
[C2(a, 6)]^, mimicking the relation between the unsubtracted functions C2 aiid C^3- Clearly 
to check for the presence of a nonzero three-point function in the COBE data, one should 
look at the smallest feasible values of a. If a signal is found, the specific prediction of 
equation (14) could be checked by looking at other configurations of a, b, and c, such as 
all three in a line on the sky. 

The nonzero three-point function implies that the fiuctuations are nongaussian. The 
particular form of the three-point function specified by equation (9) also arises (to leading 
order in A) in models in which fiuctuations in the gravitational potential (p are given by a 
local, nonlinear transformation of a gaussian random field ^: (pi^x) = ^{x) + X$,{x)'^ + . . . . 
Such models have been investigated by Kofman et al. (1991), Moscardini et al. (1991), and 
Scherrer (1992). In particular, Kofman et al. (1991) have noted that such effects can arise 
due to nonlincarities in the classical evolution equations for certain inflationary models, 
even if the initial fluctuations are gaussian. Other nongaussian models (which do not result 
from inflation) give different three-point functions. For example, Weinberg & Cole (1992) 
have studied models in which the density fluctuations Sp/p are given by a local, nonlinear 
transformation of a gaussian random fleld; this yields F oc (A:i/c2^3)^(^jf ^ + ^2^^ + ^3^^) 
to leading order in A. Still another type of model has been proposed by Scherrer & 
Schaefer (1992); it gives F oc (/ci/c2^3)'^^^- These results for F imply that, if we normalize 
the deviations from gaussian statistics at small scales, inflationary models will produce 
larger deviations on large scales than these alternatives. We would argue that future work 
on nongaussian models should focus on those which reproduce the three-point function 
reported here, since it is a generic prediction of inflationary cosmologies. 

We would like to thank Mark Wise for pointing out to us both the feasibility and 
potential interest of this calculation, and for providing us a copy of his unpublished notes 
on the related calculations reported by Allen et al. This work was supported in part by 
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FIGURE CAPTIONS 



Fig. 1. — The dimensionless two-point correlation C2(a,6) as a function of a, where a-b = 
cos a, for a = 3.2°. {Description: C2(a) is the same as the center of the gray band in 
Fig. 3 of Smoot et al. but with a different vertical scale; C2(0) = 3.04.) 

Fig. 2. — The dimensionless three-point correlation C^{a,b,c) as a function of a, where 
a-b = b-c = c-a = cos a, for a = 3.2°. {Description: C^{a) is defined only for a between 
0° and 120°. Over this range it looks a lot like C2{a)'^; Cs{0) = 9.15 and ^3(120°) = 0.) 

Fig. 3. — Same as Figure 2, with the vertical scale expanded to show the large angle struc- 
ture. 
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